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ABSTRACT 

We calculate the flux and spectrum of 7-rays emitted by a two-temperature 
advection-dominated accretion flow (ADAF) around a black hole. The 7- 
rays are from the decay of neutral pions produced through proton-proton 
^ ■ collisions. We discuss both thermal and power-law distributions of proton 

^ ■ energies and show that the 7-ray spectra in the two cases are very different. 

O ■ We apply the calculations to the 7-ray source, 2EG J1746-2852, detected by 

^ ■ EGRET from the direction of the Galactic Center. We show that the flux and 

■ spectrum of this source are consistent with emission from an ADAF around 

the supermassive accreting black hole Sgr A* if the proton distribution is 
Q_|. a power-law. The model uses accretion parameters within the range made 

Q . likely by other considerations. If this model is correct, it provides evidence 

for the presence of a two temperature plasma in Sgr A*, and predicts 7-ray 
. fluxes from other accreting black holes which could be observed with more 

sensitive detectors. 

1. Introduction 



Advection-dominated accretion flows (ADAFs) are optically thin hot accretion flows 
with low radiative efficiency (Narayan & Yi 1994, 1995 a,b; Abramowicz et al. 1995). 
Unlike standard thin disks (see Frank et al. 1992) in which the viscously generated 
energy is thermalized and radiated locally, ADAFs store most of the viscous energy and 
advect it into the central star. The gas in ADAFs has a two temperature structure 
(Shapiro, Lightman, & Eardley 1976; Rees et al. 1982), with the ions being hotter than 
the electrons. The viscous heating affects mainly the ions, the more massive species, 
while the radiation is produced primarily by the electrons. Since the ions transfer only a 
small fraction of their energy to the electrons via Coulomb scattering, the energy which 
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is radiated is much less than the total energy released during accretion (Rees et al. 
1982). 

The emission spectrum of an ADAF is mainly determined by the cooling processes 
of the electrons, viz. synchrotron, bremsstrahlung, and Compton processes (see eg. 
Narayan & Yi 1995b; Mahadevan 1997). Detailed calculations (e.g. Narayan 1996) 
show that ADAFs around accreting black holes have a characteristic spectrum ranging 
from radio frequencies ~ 10^'^ Hz up to hard X-ray frequencies ~ 10^^ Hz, whose shape 
is a function primarily of the mass accretion rate, and to some extent the mass of the 
accreting star. A number of accreting black hole systems have been shown to contain 
ADAFs (e.g. Narayan, Yi, & Mahadevan 1995; Narayan, McChntock, & Yi 1996; Lasota 
et al. 1996, Fabian & Rees 1995, Mahadevan 1997, Narayan, Barret, & McClintock 1997, 
Reynolds et al. 1997). 

In the work done so far, on ADAFs, only the cooling of the electrons has been con- 
sidered in calculating the spectra. However, since ADAFs are two temperature plasmas 
with very high ion temperatures, e.g. Tj ~ 10^^ K close to a black hole (Narayan & 
Yi 1995b), one wonders if there might be radiative processes associated directly with 
the ions. In particular, at such high temperatures, collisions between protons can lead 
to substantial production of neutral pions, 7r°, which would decay into high energy 
~ 70MeV 7-rays. Charged pions can also be produced by proton-proton collisions, but 
the final decay products are neutrinos and electrons (and their anti-particles), so that 
most of the energy leaves the plasma either in neutrinos or lower energy photons. 

Gamma-ray emission from a two temperature gas accreting onto a black hole 
has been computed previously (see e.g. Dahlbacka, Chapline, & Weaver 1974; Colpi, 
Maraschi, & Treves 1986, Berezinsky & Dokuchaev 1990). Most of the previous work 
considered a thermal distribution of protons, and the predicted 7-ray fluxes depend sen- 
sitively on the particular accretion scenario being considered. In this paper we consider 
both thermal and non-thermal distributions of protons and focus specifically on the 
ADAF paradigm. Using the methods described by Dermer (1986ab) we calculate the 
flux and spectrum of 7-ray emission from ADAFs and show that these flows produce 
interesting levels of high energy 7-rays. In the case of the source Sagittarius A* (Sgr 
A*) at the Galactic Center, we predict a 7-ray flux which is above the sensitivity limit 
of the EGRET detector on the Compton Gamma-Ray Observatory, and we compare 
the prediction with a possible detection of this source. 

The outline of the paper is as follows. In §2 we present the reaction rate for pion 
production and describe the structure and basic equations of ADAFs. In §3 we calculate 
the 7-ray spectra from thermal and non-thermal distributions of protons in ADAFs, 
and in §4 we apply these results to EGRET observations of the Galactic Center. In §5 
we discuss other applications of the results. 
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2. General Relations 



2.1. Pion Production 



The production of 7-rays from p — p colhsions is a two step process, and has been 
worked out in detail by Stecker (1971, see also Stephens & Badhwar 1981; Dernier 
1986ab). The colliding protons produce an intermediate particle, a neutral pion, 7r°, 
which then decays into two 7-ray photons. The reaction is 

p + p^p + p + n^, (1) 

7r°^7i + 72, (2) 

where the photons, 71, 72, in general have different energies in the gas frame. The 
number of neutral pions produced is calculated by considering two protons with four 
momenta pi and p2, moving towards each other with relative velocity frei. The number 
of collisions that occur in a volume dV, for a time dt, is a frame invariant quantity, which 
in an arbitrary reference frame can be written as (see Landau & Lifshitz 1975, §12) 

dRi2 = a v^ei ^ '^^ nin2 dV dt, (3) 

hii h,2 



= canin2\J{(3i-(52f-{(3iyi(32fdVdt. (4) 
Here, ni, ^2 are the number densities of the protons, E^, E2 are their energies, a is the 

— * — # 

total cross-section for pion production, and f32 are the respective velocity parameters 
of the two particles. For a distribution of particle velocities, the total number of pions 
produced per second is just the integral of dRi2 over volume and velocities, which in 
spherical co-ordinates is 



//"y /-oo rl 

R^dR j d^i J ^72 J dcose 



X a(7i, 72, cos^) n{R, 71) n{R, 72) - 0? - 0i x 0? s-\ (5) 

where 71, 72 are the respective Lorentz factors of the two protons, cos6' = /3i ■/32/|/3i| |/32|, 
and R is the radius with respect to the accreting star. The total number of 7-rays 
produced is 2N{Et^). 

To obtain the 7-ray spectrum, wc first calculate the spectrum of pion energies, 
F[Et^) dEjr, in the frame of the observer. The observed 7-ray spectrum is then obtained 
through the relation (Stecker 1971) 

F(E^) ^2 r dE^ ^ photons s'^ GeV-^ (6) 
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where = 0.135 GeV/c^ is the mass of the pion, and -E'Trmin is the minimum pion 
energy required to produce a 7-ray with energy E^, 

2 

E . = E + ^ f7) 

The pion spectrum, F{E.„) is obtained by substituting the differential cross-section 
(icr(i?^; 7i, 72, cos6')/c/_E'^ for (7(71, 72, cos6') in equation (j^). Two models are generally 
used to determine the differential cross-section for pion production: the isobar model 
and the scaling model. In the isobar model, pions are produced by a two step process. 
Colliding protons produce an isobar, with rest mass 1.238 GeV, which subsequently de- 
cays into a proton and a pion (Lindenbaum & Sternheimer 1957; Stecker 1971). Dermer 
(1986b) has shown that the isobar model agrees with the experimental data quite well 
for proton energies near threshold: Ep ^ 4 GeV. On the other hand, for proton energies 
Ep ^ 8 GeV, the scaling model of Stephens and Badhwar (1981) represents the experi- 
mental data better (Dermer 1986b). In our calculations, we therefore use the two models 
in their respective ranges of validity and interpolate smoothly (cubic interpolation) in 
the intermediate regime. 



2.2. ADAF Equations 



To evaluate equations (^) and (^ we need the number density of protons as a 
function of energy in the accretion flow, np(i?, 7). We write np(i?, 7) as the product of 
a normalized velocity distribution, which depends on the temperature at a given radius, 
and the total number density of protons, 

np{R,^)=np{R)n4j,9p{R)], (8) 

with 

n,[^,ep{R)]d^ = l, {^-l)n,[^,ep{R)]d^ = -ep{r). (9) 

Here, 6p{R) is the dimensionless ion temperature, 6p = kTi/mpC^, at radius R. 

An ADAF is characterized by four parameters: the viscosity parameter, a (Shakura 
& Sunyaev 1973), the ratio of gas pressure to total pressure, P^dv, the mass of the central 
black hole, M, and the mass accretion rate, M. Given these parameters, the global 
structure and dynamics of ADAFs can be calculated (Narayan, Kato, & Honma 1997; 
Chen, Abramowicz, & Lasota 1997). Generally, ADAFs appear to have large values of 
a ^ 0.1 (Narayan 1996), and the uncertainty in the value is a factor of a few. Many 
published ADAF models in the literature use a = 0.3 (e.g. Narayan et al. 1996, 1997). 
If there is equipartition between magnetic and gas pressure, we expect /3adv = 0.5; this 
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is the value we favor. However, for completeness we also consider the more extreme case 
Padv = 0.95, where the magnetic pressure contributes negligibly (~ 1/20) to the total 
pressure. 

The radial structure of ADAFs is well approximated by a series of concentric spher- 
ical shells, with the properties of the gas varying as a function of radius (Narayan & Yi 
1995a). The continuity equation gives M = in R"^ v{R) p{R), where v{R) is the radial 
velocity of the gas, and p{R) is the mass density. Writing the radius in units of the 
Schwarzschild radius, r = R/Rs, where Rs = 2.95 x 10^ m cm and m = M/Mq is the 
mass of the central star in solar mass units, and assuming that the mass fraction of 
hydrogen is X = 0.75, we have 



1.9 X 10^^ m'^mr'^ 



v(r] 



= m ^mnir) 



cm 



(10) 



18 



Here m = M/M^dd is the accretion rate in Eddington units, with MEdd = 1-39 x 10 
g s~^, and n{r) is defined so as to scale out the dependence of np{r) on m and m. Given 
a, Padv, fn and m, equation (|iy) allows us to calculate np(r) provided we have v{r)/c. 
We obtain v{r)/c from the global ADAF solution of Narayan et al. (1997). 

The ion temperature Tj(r) is obtained from the total gas pressure Pg (Narayan & 
Yi 1995b), 

P.-M)4{r) = ^-^^^^, (11) 

Pi nip pe nip 

where Cs(r) is the isothermal sound speed, and pi = 1.23, pe = 1.14 are the effective 
molecular weights of the ions and electrons respectively. Since T, ^ Tg we neglect the 
second term to obtain 

2 



Ti(r) = 1.34X 10"/?adv 



K, 



(12) 



where c is the velocity of light. We obtain Cs(r)/c again from the global solution of the 
ADAF (Narayan et al. 1997). 



Given np{r) and Tj(r) and an assumed functional form for the proton energy dis- 
tribution ny{'y,9p), we can substitute equation 
the 7-ray spectrum. 



in equations (^ and @ to calculate 



3. Results 



At present our understanding of viscous heating is poor, and we do not know 
whether the process leads to a thermal or power-law distribution of proton energies, or 
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perhaps some combination of the two. In this section we calculate the pion spectra by 
considering two different proton energy distributions: a relativistic Maxwell-Boltzmann 
distribution and a power-law distribution. 



3.1. Thermal Distribution 



The proton temperature in an ADAF close to a black hole is marginally relativistic, 
9p = kTp/nip ^ 0.2. At such temperatures, pions are produced primarily by protons in 
the tail of the Maxwell-Boltzmann distribution, since only these particles have energies 
above the threshold needed for pion production. Thus, the pion production is very 
sensitive to the proton temperature, and since 6p scales approximately as 1/r (Narayan 
& Yi 1995b; Mahadevan 1997), most of the production is from the innermost radii of 
the ADAF, r < 10. 

The protons at these temperatures have momenta near the threshold for pion pro- 
duction, and do not exceed ~ 1 GeV/c. At such energies the isobar model (Lindenbaum 
& Sternheimer 1957; Stecker 1971; Dermer 1986a) is more accurate and we have therefore 
used it to calculate the cross-sections and energy spectra. In this model, the isobars 
move along the initial directions of the colliding protons in the center of momentum 
(CM) system. It is therefore convenient to transform variables in equation (|^) via a 
Lorentz boost to the CM reference frame, followed by a rotation so that the colliding 
protons are moving along the z-axis (Dermer 1984, 1986a). Using equation (^ with 
equation ([10|) , and rewriting in terms of the new variables (see Dermer 1984, 1986a for 
details), we obtain 

(2.95 X 10^)^ mm^ TTc T™^" , r^n^(r) i"^ (7^-1) 



'2 



Ji ep{r)Kl[l/ep{r)] Ji [2(7, + l)]i/ 

X Z^™^:^^^^^^ {exp[-g77*(l-/?/5*)]-exp[-g77*(l + /?/?*)]}, 

photons s-i GeV-^ (13) 

Here E^^ = rriT^'j is the pion energy in the observer's frame, 7^ = 7172(1 — /?i/?2 cos6') 
is the relative Lorentz factor of the two colliding protons, q = [2(7r + l)]^^^/^p(r), 
da*{'y*; 'yr)/d'y* is the differential cross-section for the production of a pion with Lorentz 
factor 7* in the CM frame for a given 7^, C = {S — Arrip+m'^) / S^^'^ where S = 2771^(7^ + 1) 
characterizes the strength of the collision, K2{l/6p) is the modified Bessel function of 
order 2, and we have used a Maxwell-Boltzmann energy distribution for the protons. We 
take Tinax to be 10^; however, since the proton temperature decreases steeply with radius 
{Ti ~ 1/r) the contribution to the pion spectrum from the gas at r > 30 is negligible. 



The calculation in equation (|T3D gives the pion spectrum, F{Et,). Using equation 



we then obtain the 7-ray spectrum. The results are shown in Fig. 1. The plot 
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corresponds to two values of the viscosity parameter a, 0.1 and 0.3, and two values 
of /5adv5 0.5 and 0.95. For a given a, increasing /5adv leads to increasing gas pressure, 
and therefore to increasing ion temperature (see eq. O). Changing /5adv from 0.5 to 
0.95 causes the temperature to go up by a factor ~ 2, and since the pion production is 
extremely sensitive to temperature, the 7-ray luminosity increases by nearly two orders 
in magnitude. Changes in a also affect the 7-ray luminosity, though less sensitively. 
Since the total luminosity is proportional to r?^ which varies roughly as (Narayan 
& Yi 1995b), we might expect the luminosity to increase with decreasing a. This effect 
is counterbalanced however by small changes in the temperature, and the luminosity in 
fact decreases slightly when a goes down from 0.3 to 0.1. 



3.2. Non— Thermal Distribution 



In this section we consider a power-law distribution of proton energies described 
by a spectral index s. Equation (|l^) continues to hold, but n[7,^p(r)] is no longer a 
Maxwell-Boltzmann distribution. To determine the form of n[7, ^p(r)], two requirements 
have to be satisfied (see eq. P), namely the normalization condition and the requirement 
that the average kinetic energy of the power-law distribution of protons, at each radius, 
be equal to the average energy of the protons as given by the local ion temperature. 

We model the energy distribution of protons as 

n[7, e,{T)\ ^7 = { [1 - C(r)] 5(7 - 1) + (s - 1) C(r) 7"^} ^^7, (14) 

where a fraction 1 — C(r) of the protons have 7 ~ 1, and a fraction C^ir) are in a power- 
law tail with index s. The distribution is properly normalized. The fraction C(r) is fixed 
by the energy requirement: 

C(r) = ^(s-2)e,(r). (15) 

Since the exact energy distribution for the fraction \ — (^{r) of protons with energy below 
threshold is not important for the present calculation, we have simplified the distribution 
to a (5 function in equation (|T3|). 

To calculate the total pion luminosity, N^E^,), we substitute equation ( p^ in equa- 
tion (^. We then obtain three different terms proportional to (1 — C)^, (1 — () (, and 
C^, respectively. The (1 — C)^ term gives no contribution since it corresponds to protons 
with 7 ~ 1 colliding with one another. These protons do not have enough energy to 
produce pions. The (1 — C) C term corresponds to protons with 7 > 1 colliding with 
protons with 7 ~ 1, and we denote the contribution by NiIEt^). Similarly, the term 
corresponds to protons with 7 > 1 colliding with protons with 7 > 1, and we denote it 
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as N2{E^). We have 

Ni{E^) = 4:nc{2.95xl0ymm'^{s-l)Jdrr^h\r)[l-C{r)]C{r) 

/oo 
CTmhilh'' [3 d-f s-\ (16) 

and 

drr^n^{r)C\r) J dji-f^'J ^72 72"' 
rf(cos^) cr^b(7i, 72, cos^) \/ {Pi - /S^Y - {Pi X s-\ (17) 



1-1 

where cxmb is the cross-section in units of milhbarns, and 

„2 -2/ 



J(s,a,/3adv) = 9.68(s-l) J drr'n'{r)[l-C{r)]C{r). (18) 

In an ADAF, we generally have ( <^ 1 so that most of the protons have 7 ~ 1 (cf. eq. 
p!5|). We therefore expect most of the pion production and 7-ray luminosity to come 
from Ni{E.,,). Table |l| gives numerical results and we see that N2{Et,) is indeed negligible 
compared with Ni{Et,). 

In calculating the pion spectrum, we consider only the contribution from Ni{Et^) 
and write 



dE^ 

J(s,a,/3adv)mm2 _5 ^ cio-jnb(7,r,7) 



r d^r'P photons s-i GeV-i.(19) 

Jl d'y^ 



rriT^ Jl d'jTj 

Unlike in the case of a thermal distribution where most of the protons have energies 
either below or just above threshold, in a power-law distribution the protons have a 
wide range of energy and some protons are well above threshold. We therefore use both 



the isobar and scaling models, as described in § |2.1 



Equation ([19D reveals the dependence of the pion flux and energy spectrum on the 
parameters: (1) The flux is proportional to mrn?. (2) The flux depends on a, /5adv 
and s through the function /(s,a,/3adv) defined in equation (|T8|). (3) The shape of the 
spectrum depends only on s and is given by the integral in equation ([I9|). 

Tables ^ and ^ give /(s, a, /5adv) for different values of its parameters. For fixed 
a, increasing /5adv leads to increasing /(s,a,/5adv) since the temperature in the flow 
increases. However, unlike the thermal case, the luminosity is not excessively sensitive 
to /3adv since the number of protons above threshold is directly proportional to 6p (cf. 
eq. |1^); therefore, a change in /3adv by a factor of two changes the total luminosity only 
by ~ 2. For fixed /?adv in a self-similar solution (Narayan & Yi 1994, 1995b), the density 
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varies in proportion to a~^, so the luminosity scales as However, in the more 

accurate global solutions employed in this paper, the dependence of density — and hence 
luminosity — on a is not strictly a power-law, and the overall strength of the dependence 
is somewhat weaker than in the self-similar solutions. To illustrate this dependence, we 
list the numerical values of /(s,a,/3adv) for three different choices of a in Tables ^ and 
I- 

Figure 2 shows 7-ray spectra for different values of the proton energy spectral index 
s (we fix tt = 0.3, /3adv = 0.5). The 7-ray spectral index at high energies 1 GeV 

is the same as the energy spectral index of the protons s, as shown by Dermer (1986b). 
For this reason, we investigated values of s in the range 2.3 - 3.3. The spectrum turns 
over at ~ 70 MeV (approximately half the pion rest mass) and falls for lower photon 
energies. Comparing Figures 1 and 2, we see that the total 7-ray luminosity from a 
power-law energy distribution of protons is comparable to that obtained from a thermal 
distribution when /3adv = 0.95, but is very much higher for /?adv = 0.5. In our view 
Psidv = 0.5 is the natural choice since it corresponds to equipartition between gas and 
magnetic pressure. The 7-ray spectrum extends to much higher photon energies when 
the protons have a power-law distribution compared to the thermal case. 



4. Application to Sgr A* 

The black hole candidate at the center of our Galaxy, Sagittarius A* (Sgr A*), 
appears to be a scaled down version of an AGN. The source is believed to consist of a 
black hole of mass M = 2.45 x lO^M© (Eckart & Genzel 1996a, 1996b) accreting gas 
from its surroundings. 

The mass accretion rate in Sgr A* has been estimated by various methods (see 
Genzel et al. 1994 for a summary). The accretion of wind gas from neighboring stars 
(principally IRS 16) seems to be the most likely scenario. The accretion rate is deter- 
mined by considering the fraction of the wind that comes within the accretion radius 
of the central black hole, but the estimate depends sensitively on the assumed wind 
velocity. Genzel et al. (1994) obtain an accretion rate of M ^ 6 x 10~^ ^0/ yr~^ which 
corresponds to m ~ 10"'^ for M = 2.45 x 10^ Mq, assuming a wind velocity of 1000 
km s~^ and mass loss rate from winds of ~ 3.5 x 10"'^ ^0/ yr~^- Melia (1992) uses a 
wind velocity of 600 km s~^ and obtains an accretion rate of M ~ 2 x 10~^ -^0/ yr~^ 
(m ~ 3 X 10~^). The true accretion rate probably lies between these two estimates, and 
we will require our models to have 10~^ < m < 3 x 10^^. Note that Lacy et al. (1982; 
see also Rees 1982) estimated a much higher mass accretion rate of M ~ 2 x 10~^ Mq/ 
j^-i fj^Qj2i stellar disruptions at the Galactic Center, but argued that the accretion due 
to this probably has a low duty cycle. 
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Sgr A* has an extremely low luminosity, L ~ 10^^ ergs~^, which corresponds to 
m ~ 10~^ if the accretion flow has a standard radiative efficiency of 10%. Such a 
low m is in serious conflict with the estimate of fn given in the previous paragraph. 
Narayan, Yi & Mahadevan (1995) suggested that Sgr A* does accrete at a rate near 
the one estimated, but in an advection-dominated mode. Using a self-similar ADAF 
model, they were able to reconcile the low luminosity of the source with a relatively 
large m: ~ 8 x 10~^q;. If q; > 0.1, the accretion rate inferred on the basis of the 
ADAF model would be in agreement with that estimated on the basis of gas supply. 
Furthermore, they obtained a reasonable fit to the observed spectrum from radio to 
hard X-ray frequencies. In this section, we calculate the 7-ray luminosity and spectrum 
of Sgr A* due to pion production, and compare these predictions with the fiux seen by 
EGRET from the direction of the Galactic Center. 

Among several unidentified sources detected by EGRET in the Galactic plane 
(Merck et al. 1996), the source 2EG 1746-2852 is of particular interest since it is spa- 
tially coincident with the Galactic Center. There is no firm identiflcation of this source, 
but it is point-like to within the resolution of the instrument (~ 1°), it is ~ lOo" above 
the local diffuse emission, and its spectrum differs significantly from the spectra of other 
unidentified EGRET sources. Out of the 32 sources whose spectra are reported by Merck 
et al. (1996), 27 have spectral slopes > 2.0 (photon index), 4 have spectral slopes ~ 1.9 
and one has a very hard spectrum with a spectral slope of 1.7. This last source is 2EG 
1746-2852, and it is located exactly at the Galactic Center. Following Merck et al. 
(1996), we make the reasonable assumption that 2EG 1746-2852 corresponds to Sgr A*. 
Figure 3 shows the spectrum of the source as measured by EGRET. The dashed line is 
the best-fit power-law obtained by Merck et al. (1996) over the photon energy range 
100-4000 MeV. Their fitted spectral slope is s = 1.7 ± 0.1. At higher energies, the data 
suggest a roll-over in the spectrum which would give a softer spectral index at higher 
energies (E 3> 4 GeV). 



4.1. Thermal Distribution 



Figure 3 compares the 7-ray spectrum from the thermal model with the EGRET 
data. The solid fine corresponds to a = 0.3, /3adv — 0.95, and the dotted line to a = 0.3, 
/3adv = 0.5. The black hole mass is taken to be m = 2.45 x 10^, and the accretion rates 
in the two models have been varied such that the predicted spectra agree with the data 
point at £^ ~ 200MeV; the accretion rates are indicated on the plot. 

It is clear that the thermal model predicts a spectrum with the wrong shape and 
fails to explain the EGRET detection at £■ ~ 1 GeV. Moreover, for the preferred value 
of /3adv = 0.5, the model requires a rather high m, which is much larger than the range 
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we consider reasonable (see above). We conclude that the thermal model is inconsistent 
with the observations. 



4.2. Non— Thermal Distribution 

Fig. 4 compares 7-ray spectra from various power-law models (a = 0.3, /5adv = 0.5) 
with the EGRET data. The different curves correspond to different values of the proton 
energy index s. For each s, we have varied the accretion rate so as to obatin the closest 
agreement with the observed spectrum. The results (5 x 10~^ < rh < 9.2 x 10~^), which 
are only weakly dependent on s, are given in Table ^. The accretion rates so inferred 
fall well within the range of accretion rates estimated on the basis of the properties of 
nearby gas (10"'^ < m < 3 x 10^^, Genzel et al. 1994), and are within factors of a few 
of the accretion rate estimated on the basis of fitting the radio to X-ray spectrum with 
the somewhat less accurate self-similar version of the ADAF model (m ~ 8 x 10~^a, 
Narayan, Yi, & Mahadevan 1995). 

Reproducing the power-law portion of the 7-ray spectrum is a by-product of pos- 
tulating a power-law proton energy distribution. However, it is striking in this context 
that the best-fit slope, s ^ 2.7, is identical to the slope of the low-energy cosmic ray 
distribution. The models predict roll-overs at both the high energy (few GeV) and low 
energy (~ 200MeV) ends of the spectrum. At the high-energy end there is just a hint 
of such a roll-over in the data, but at the low-energy end, the evidence seems somewhat 
stronger. 

The similarity of the best-fit spectral shape to the diffuse Galactic 7-ray emission 
has two possible interpretations. First, since the diffuse emission is generally interpreted 
as being due to cosmic ray protons striking thermal protons in the interstellar medium, 
a mechanism similar to what occurs in an ADAF, the EGRET detection of 2EG 1746- 
2852 could be interpreted as cosmic ray protons colliding with a compact dense cloud 
of gas that is spatially coincident with the Galactic Center. This interpretation would 
require such a cloud to have a considerably greater gas density (and possibly also cosmic 
ray density) than in the interstellar medium only slightly farther away from the Galactic 
Center. If virtually all the Galactic Center region 7-ray flux were due to such a cloud, 
the power-law proton distribution function version of the ADAF model would be put 
seriously in doubt. Thermal proton distribution function versions of the ADAF model 
(or entirely different models) would not be seriously constrained if the 7-rays coming 
from the Galactic Center prove not to have their source in Sgr A*. 

The alternative is that the observed 7-ray spectrum is in fact due to Sgr A*. In this 
case, not only is the deduced accretion rate roughly consistent with previous estimates, 
but we also have the interesting result that whatever process determines the cosmic 
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proton energy distribution (shock acceleration ?) is also at work in ADAFs. We find 
the numerical coincidence between the prediction of the power-law version of the ADAF 
model and the observed 7-ray flux striking enough to justify pursuing this interpretation. 

To gauge just how strong the quantitative agreement is, we must consider how 
much freedom the model is given by adjustable parameters. Three free parameters are 
significant. One (m) is fixed to within an order of magnitude or so by other considerations 
(the amount of gas available to accrete, and fitting the radio-X-ray spectrum). Another 
(a), while formally unconstrained over a range of several orders of magnitude, is often 
supposed to have a value quite near the one which we infer (i.e. between 0.1 and 1, cf. 
Narayan 1996). Only the third (s) is chosen almost entirely on the basis of fitting the 
observed 7-ray spectrum. Because the 7-ray luminosity predicted by the ADAF model 
(assuming a power-law proton energy distribution) approximately scales with 'm^a~^^^ 
(the actual scaling with a is only very roughly described by a power-law) , the combined 
a priori uncertainty in m and a may be regarded as giving the 7-ray luminosity predicted 
by the model a possible range of ~ 100. 



5. Discussion &c Conclusions 

Since it is not understood whether viscous heating produces a thermal or power- 
law distribution of proton energies, we have calculated the 7-ray emission of ADAFs 
corresponding to both types of distributions. We expect that the true energy distribution 
in any given source will be bracketed by these two extremes. Spectra corresponding to 
intermediate models can be easily calculated by taking a weighted sum of the thermal 
and power-law models. 

If the proton distribution is thermal, the spectrum has a characteristic shape with a 
peak at ~ 70 MeV and very little emission at either lower or higher photon energies. 
For the kinds of proton temperatures expected in an ADAF, only a small fraction of the 
protons (in the tail of the Maxwellian distribution) have sufficient energy to produce 
pions. Consequently, even a minor change in the temperature, say by a factor of ~ 2, 
can modify the 7-ray luminosity by orders of magnitude (cf. Figure 1). Since the flux 
is very sensitive to the gas temperature, a detailed understanding of the physics of the 
ADAF in the region r ^ 10 (where most of the 7-ray luminosity originates) is necessary 
in order to make testable predictions. The present study does not have the necessary 
accuracy for this. General relativistic effects become important at these radii and need 
to be included consistently. In principle, if these effects are included, the 7-ray spectrum 
could be used as a sensitive probe of the proton temperature. One might even hope to 
distinguish between rotating and non-rotating black holes, since the physical properties 
of the flow at r ^ 10 will be different for the two cases. 
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If the protons have a power-law distribution of energies, the 7-ray spectrum again 
peaks at ~ 70 MeV, but there is significant emission at higher photon energies. Indeed, 
the spectrum asymptotically has a spectral slope s which is equal to the power-law 
index s of the proton distribution (Dermer 1986b). Thus, the detection of a power-law 
7-ray spectrum not only indicates the nonthermal nature of the protons but also helps 
determine the energy index. 

In contrast to the thermal case, only half the 7-ray luminosity due to a power-law 
distribution of proton energies originates from r ^ 10, while the other half is emitted 
from r ^ 10. Further, the luminosity is not very sensitive to changes in the temperature. 
Therefore, our present understanding of the physics of ADAFs is probably adequate for 
a reasonable estimate of the 7-ray flux, and we can thus test this version of the ADAF 
paradigm usefully against observations. 

The Galactic Center source, Sgr A* presents an excellent opportunity to test the 
model, since there is a good case for the presence of an ADAF in this source (Narayan 
et al. 1995) and the predicted 7-ray flux is large enough to be detectable with current 
instruments. Indeed, EGRET has detected an unresolved source in the Galactic Genter 
region, 2EG J1746-2852 (Merck et al. 1996). It is possible that this source corresponds 
to a dense compact gas cloud interacting with cosmic rays. However, the source could 
equally well be associated with Sgr A*. The 7-ray flux of Sgr A* which we calculate with 
our model using a power-law distribution of proton energies is in good agreement with 
the observed flux of 2EG J 1746-2852 given parameters which are almost determined 
on independent grounds. In order to fit the fiux, we need an Eddington-scaled mass 
accretion rate of m ~ 5-9 x 10~^ (for a — 0.3), close to the estimate m ~ 8 x 10~^q; 
obtained by Narayan et al. (1995) from fitting the lower energy spectrum due to the 
electrons. It is also compatible with the range of m quoted by Genzel et al. (1994) on 
the basis of gas motions in the vicinity of Sgr A*: 10"^ < m < 3 x 10~^. In addition, the 
pion-decay model for the Galactic Genter 7-ray spectrum predicts a roll over at energies 
below ~ 100 MeV that may be present in the observed spectrum. 

An exciting aspect of the present study is that it provides for the first time a direct 
probe of the protons in hot accretion fiows. A fundamental assumption in ADAF models 
is that the plasma is two-temperature. However, until now there has been no direct test 
of this assumption since all previous investigations of the emission from ADAFs dealt 
only with the electrons. The 7-ray emission we have considered in this paper is due 
entirely to the protons. Moreover, it requires that the protons have nearly virial energies 
in order to be able to exceed the energy threshold for pion production. The encouraging 
results we have obtained in the case of Sgr A* indicate that this source might have a 
two-temperature plasma exactly as postulated in the models. More detailed study of 
Sgr A* coupled with future detections of other sources (see below), could strengthen the 
case for two-temperature accretion fiows significantly. 

Another exciting aspect of this study is that the 7-ray spectrum from pion de- 
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cay provides a direct window to the energy distribution of the protons. In the ADAF 
paradigm, viscous heat energy goes almost entirely into the protons. Furthermore, it is 
easy to show that, at least at the low m expected in quiescent systems, the protons have 
no energetically important interactions either among themselves or with the electrons. 
Therefore, each proton retains memory of all the heating events it has undergone during 
the accretion, and so the energy distribution of the protons directly reflects the heating 
processes present in the plasma. In principle, with sufficiently sensitive observations, 
one could use 7-ray measurements as a direct probe of viscous heating. This would be 
invaluable for the theory of hot accretion flows. 

The ADAF model has been applied successfully to several other low-luminosity 
black holes in addition to Sgr A*. The sources studied include two X-ray binaries in 
quiescence, A0620-00 and V404 Cyg (Narayan, McClintock, & Yi 1996, Narayan, Barret 
& McClintock 1997), and several quiescent AGNs, viz. the LINER galaxy NGC 4258 
(Lasota et al. 1996, Herrnstein et al. 1996) and the nearby elliptical galaxies NGC 
4472, NGC 4486, NGC 4649, and NGC 4636 (Mahadevan 1997; Reynolds et al. 1997). 
Table ^ gives the values of m, m, a and /3adv of several of these systems (taken from 
the references indicated) and presents the expected 7-ray fluxes according to our model. 
Among the nearby ellipticals, we have included only NGC4486 since it is the only system 
with a reliable mass estimate (Ford et al. 1995; Harms et al. 1995). The 7-ray flux 
estimates in Table |^ assume that the protons have a power-law energy distribution with 
s = 2.75. We see that, only in the case of Sgr A* does the detection threshold of EGRET 
(10"*^ photons cm~^ s^^) permit a test of the predictions, although some of the other 
sources might be not far below the EGRET threshold. These sources are potentially 
detectable with future instruments such as the Gamma Ray Large Area Space Telescope 
(GLAST), which is designed to be 100 times more sensitive than EGRET. 

Narayan (1996) has argued that X-ray binaries in the "hard state" (also called the 
"low state") may also contain ADAFs with rh ^ 0.1. If the proton energy distribution 
in these sources is a power-law of the sort we suggest exists in Sgr A*, several of them 
might be expected to produce detectable 7-ray fluxes. However, the accretion rates 
in these objects (e.g. Cyg X-1) is likely to be close to the maximum rate permitting 
an ADAF (Narayan & Yi 1995b). If this is the case, particle interactions might be 
rapid enough that maintenance of a true power-law distribution function would be 
questionable. Accurate predictions of the 7-ray flux would then require a more elaborate 
calculation than we have performed for Sgr A*. 

Finally, we note that in determining the 7-ray spectra, we have neglected gravita- 
tional redshift effects which become important at r ^ 5. We have also neglected the 
Doppler blueshift associated with the large radial and orbital velocities of the gas at 
these radii (Narayan & Yi 1995b). Only a detailed calculation can tell which of the two 
effects predominates. We note, however, that these effects are not likely to modify the 
results we have presented for a power-law distribution of protons because more than 
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half the emission in this case occurs at r ^ 10. 
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Table 1: Photon luminosities (photons s ^) in units of mrn^ for a = 0.3, /3adv = 0.5, and 
different power-law indices s. 



s 




N2{E^) 


N{E,) 




mm? 


mrn? 


mrn? 


2.1 


7.36e+39 


1.52e+37 


7.38e+39 


2.2 


1.26e+40 


5.27e+37 


1.27e+40 


2.3 


1.65e+40 


1.04e+38 


1.66e+40 


2.4 


1.92e+40 


1.62e+38 


1.94e+40 


2.5 


2.12e+40 


2.23e+38 


2.14e+40 


2.6 


2.26C+40 


2.86e+38 


2.29e+40 


2.7 


2.36e+40 


3.48e+38 


2.39e+40 


2.8 


2.42e+40 


4.11e+38 


2.46e+40 


2.9 


2.47e+40 


4.61e+38 


2.55e+40 


3.0 


2.49e+40 


5.23e+38 


2.54e+40 


3.1 


2.50e+40 


5.72e+38 


2.55e+40 


3.2 


2.49e+40 


6.15e+38 


2.55e+40 


3.3 


2.47C+40 


6.80e+38 


2.54C+40 


3.4 


2.45e+40 


7.01e+38 


2.52C+40 


3.5 


2.42e+40 


7.50e+38 


2.49e+40 
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Table 2: Selected values of I{s, a, /3adv)/10^° for /3adv = 0.5 



a 


s = 2.1 


s = 2.3 


s = 2.75 


s = 3.3 


0.03 


0.46 


1.6 


5.4 


12. 


0.1 


.11 


0.40 


1.3 


3.0 


0.3 


0.024 


0.082 


0.28 


0.62 



Table 3: Selected values of I{s, a, /?adv)/10^° for /?adv — 0.95 



a 


s = 2.1 


s = 2.3 


s = 2.75 


s = 3.3 


0.03 


0.88 


3.2 


10. 


24 


0.1 


.15 


.52 


1.7 


3.8 


0.3 


0.036 


0.13 


0.42 


0.94 



Table 4: Mass accretion rates required in Sgr A* in order to fit the observed 7-ray 
spectrum. The models assume m — 2.45 x 10^, a — 0.3,/3adv = 0.5, and a power-law 
distribution of protons with index s. 



s m{m = 2.45 x 10^) M (Mg/yr) 

2I 9.2 xlO-^ 5.2 xlO-^ 

2.3 6.6 xlO-^ 3.5 xlO"^ 

2.75 5.2 xlO"^ 2.8 xlO'^ 

3.3 4.7 xlO"^ 2.5 xlO'^ 
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Table 5: Flux of photons > 100 MeV from various accreting black holes with ADAFs (s 
= 2.75). 



Name 












Flux (photons cm ^ s ^) 


a 




m 


m 


-Dkpc 


Sgr A* 


0.3 


0.5 


2.45 xlO^ 


5.2 xlO"^ 


8.5 


4.9 xlO"^ 


A0620-00 


0.3 


0.5 


6 


1.2x10-3 


1 


2.8x10-*^ 


V404Cyg 


0.3 


0.5 


12 


4.6x10-=^ 


3 


1.5x10-*^ 


NGC 4486 


0.3 


0.5 


3x10^ 


10-2.5 


16x10^ 


6x10-9 


NGC 1258 


0.3 


0.95 


3.6x10' 


10-- 


6.5x10-^ 


LlxlO-"" 
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Figure Captions. 

Figure 1: 7-ray spectra from a thermal distribution of colliding protons for two different 
values each of a and /3adv 

Figure 2: 7-ray spectra from a power-law distribution of colliding protons for different 
values of the energy index s. The curves correspond to a = 0.3, /5adv = 0.5. Tables || 
and 1^ can be used to obtain fluxes for other values of a and /3adv 

Figure 3: 7-ray spectra from a thermal distribution of colliding protons in the Galactic 
Center source Sgr A*. The data correspond to EGRET observations of the unidentified 
source 2EG J1746-2852, which has been identified with Sgr A* (Merck et al. 1996). 

Figure 4: Similar to Figure 3, but for a power-law distributions of proton energy. The 
curves correspond to different values of the proton spectral index s = 2.1, 2.3, and 2.75. 
Note the significantly superior fit compared to Figure 3. The mass accretion rates of the 
various models are listed in Table 4. 



Figure 1 



40 



1 1 1 — 

a = 0.3 
a = 0.1 



1 1 r 



"1 1 r 



I 

0) 38 



Pi 

O 



36 



34 



32 








2 



4 



Log(E^ MeV) 



6 



Figure 2 




2 4 

Log(E^ MeV) 



Figure 3 



-6 



8 



> 




m 
I 



o 

in^ 1 
O 

-1-3 

o 

X 

1-12 

O 



-14 







1 1 1 1 1 1 1 1 1 1 r 

a = 0.3, (3^^^ - 0.95, m = 7.2 xlO-^ 
a = 0.3, = 0-5, m = 8.0 xlO-3 




2EG J1746-2852 

EGRET 



2 



4 



6 



Log(E^ MeV) 



Figure 4 




